Lattice measurement of Bb s with a chiral light quark action 
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The computation on the lattice of the bag parameter Bb b associated to the B s — B a mixing amplitude in the 
Standard Model is presented. The estimation has been made by combining the static limit of HQET and the 
Neuberger light quark action which preserves the chiral symmetry on the lattice. We find Bj s Sstat (m(,) = 0.92(3). 



1. Introduction 

B s — B a mixing is highly important in testing 
the Standard Model (SM) and constrains strongly 
its extensions by bounding the unitarity triangle. 
Since it is a flavor changing neutral process it oc- 
curs through loops so that the corresponding mix- 
ing amplitude is a sensitive measure of \Vt s \ and 
\Vtb\, as the major SM loop contribution comes 
from t-quark. The mixing of weak interaction 
eigenstates B s and B s induces a mass gap AM, 
between the mass eigenstates B s h and B s l- Ex- 
perimentally, DO has bounded AM S as 17 ps _1 < 
AM S < 21 ps" 1 (90% CL) [I] and the measure- 
ment made at CDF gave AM S = 17.3301^ 0. 
Theoretically the B s — B s mixing is described by 
means of an Operator Product Expansion, i.e. the 
Standard Model Lagrangian Csm is reduced to 
an effective Hamiltonian H^ =2 , up to negligible 
terms of 0(1 /M^): 

nti= 2 = ^M 2 w (V t tV ts ) 2 

x r, B s (x t )C(n b )QtE =2 (^), (i) 

where t\b = 0.55 ± 0.01. So(xt) is a known 
Inami-Lim function of x t = m^/M^ [3], C{n b ) 
is the Wilson coefficient computed perturbativcly 
to NLO in a s (fj,b) in the MS (NDR) scheme, and 
Qll~ 2 IS a four-fcrmions operator coming from 
the reduction of the box diagrams in Csm to a 
local operator in the effective theory. The mass 
splitting is 

AM B , = 2\(B s \H^ 2 \B)\ . (2) 



The hadronic matrix element of QlE =2 , which 
must be computed non perturbatively, is conven- 
tionally parameterized as 

{Bl\Q$l= 2 {pi b )\B s ) ee \m B J 2 B B B M , (3) 

where _Be s (/Xb) is the B s meson bag parameter 
and /b s the decay constant. In the following it 
will be useful to introduce in addition to the op- 
erator Oi = Qll =2 the operators of the so called 
supcrsymmetric basis 

2 = 6 t (l-7 5 )s l 6 J (l- 7 5 )s J , 

3 - 6^(l-7 5 )^7M(l+7 5 y, 

Oi = ^(l-7 5 )s ?; 6 J '(l+7 5 )s J , (4) 

whose the matrix elements (B s \Oi\B s ) are param- 
eterised by 

{Bs\0 2 \B s ) = -\( V f 2 B BM 

(Bs\0 3 \B s ) = -i ( ) f B B 3 (fx), 

(b;\0 4 \B s ) = 2( ) f B B^).(5) 

So far BB 3 (fJ- b ) has been computed by using lat- 
tice qcd u-nnS One of the major problems 

with those computations is in the following: the 
standard Wilson light quark lattice action breaks 
cxplicitely the chiral symmetry, which tremen- 
dously complicates the renormalization procedure 
of Qll =2 and its matching to the continuum. To 

1 It was also estimated by using QCD sum rules but 
we will concentrate only on lattice results. 
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get around that problem we compute Bp s (//&) by 
using the lattice formulation of QCD in which the 
chiral symmetry is preserved at finite lattice spac- 
ing |12j . On the other hand, it should be stressed 
that our heavy quark is static, as the currently 
available lattices do not allow to work directly 
with the propagating b quark. Thus our results 
will suffer from l/mb-corrections. 



fi-i»0 ^ 



Z Av {B s \O x {n)\B s ) v e-'^, (7) 



C^ s+PP (t u t) = (^^+(x,i l )O 2 (0,0)A- t (y,t)) f 



2. Computation on the lattice 

In our numerical simulation we choose to work 

where 



with the action S = 



Z A v (B s \0 2 (ri\B s ) v e-'to-Q, (8) 



SF* = 



a 3 [h + (x) - V* YP \x - 6)h+(x - 0)] Q 2 = P< P^U^ P L j 



'Z A = (0|A±|B s 



h~(x) [V^ YP (x)h-(x + 6) - h~(x)] } 



is the static limit of HQET action [T3] which has 
been modified after using the so-called HYP (hy- 
percubic blocking) procedure [14] . that is enough 
to substantially improve the signal/noise ratio 
[TS] [the field h + (h~) annihilates the static heavy 



quark (antiquark)]. 



(m ) 



is the overlap light quark action with 
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where D\y is the standard Wilson-Dirac oper- 
ator. The overlap Dirac operator Z?^™ ' veri- 
fies the Ginsparg- Wilson relation {7 s , D^™ '} = 



a n W 5 n(mo) 

~ U AT J Up 



and the overlap action is invari- 



p-'N 1 

ant under the chiral light quark transformation 



ip — » ip + z£7 



, a n (mo) 

1 -o N 



which is essential to prevent mixing of four- 
fermion operators of different chirality [17] . In 
other words, in the renormalization procedure, 
the subtraction of the spurious mixing with d = 6 
operators will not be needed. We thus compute 
the two- and three-point functions: 



cf A {t) = (£,At(2,t)jtf(p))u ^ Z A e^\ (6) 



e is the binding energy of the pseudoscalar heavy- 
light meson. In the computation of C^'^{ti,t) 
one current is local whereas the other is 
smeared. The role of the smearing is to isolate 
earlier the ground state [18], as shown in Fig. 
[T][E We see that the same state is isolated when 
purely local currents are used (with those currents 
the signal does not exist if Vq HYP is not used in 
the heavy quark action). The source operators in 



*(3) 



(ti,t) and Cgg +PP (ti,t) are the smeared 



VV+AA 



a y 

currents ^4^, whereas the four-fermion operators 
Ox and O2 are purely local. In ((6]), ([7]) and ((HJ) the 
subscript " v" and superscript " ~" are designed to 
remind the reader that states and operators are 
defined in HQET. 

Note that in the computation of C vv , AA (ti,t) 

and C^g +PP (ti,t) there are two terms, coming 
from two different Wick contractions: 

^/Un^T and J^B^B^U). 

i 3 i,j 

i, j are the color indices and 



By(i)=Tr 



7lML Sl ik (0;x,th 0l 5 S k J(x,f,0) 



Sl and Sh are the light and heavy propagators 
respectively and the trace is over spinor indices. 



2 Even if the time interval from which we extract the bind- 
ing energy starts at t = 9 (green line), the overlap with 
radial excitations is quite reduced since t = 6 when cur- 
rents are smeared. 
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Figure 1. Effective binding energy of the CP- 
state when currents are local (unfilled symbols) 
or smeared (filled symbols) 
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• R,(t,=8, t) 
□ R 2 (t.=8, t) 

_u o o U $ $ 



[}=S.O, m =0.06 GeV, p=1.4, V=16 "32 



Figure 2. Signals for i?i,2 defined in eq 
lines indicate the time interval on which we fit the 
signal to a constant to extract 8i(a) and .82(a) 
respectively. 



After having computed the correlation functions 
©, and ijHJ) we build the following two ratios 
Jfcfa.t) and R 2 (ti,t): 



Ri(U,t) 



C 



(3) 

VV+AA 



{tut) 



-t»0 



„lB B |Oi|ff a 
|l(0|iol^) 



^2 (tj) i) 



C 



(3) 

SS+PP 



- = BUa), 
(*<,*) 



~\^A^AA + (^)^AA (*) 



ti-t>0 



= B 2 (a).(9) 



Those ratios are calculated either with a fixed 
time f € [-6, -8, -fO, -12, -14, -16] and t t free, 
or by fixing U € [6,8, 10, 12, 14, 16] while letting 
t free. We take the average of the two options. 



In Fig. [2] we show the quality of the signals for 
B-i 2(ti, t), with ti — 6 fixed. The signal for Bi(a) 
is quite stable as a function of ti, whereas the 
signal for 82(a) rapidly deteriorates for larger ti, 
and is completely lost for ti > 10. 



3. Extraction of physical Bb s 



Three steps are required to extract Bb b 
from the lattice: 

(1) 81.2(a) are matched onto the continuum 
MS(NDR) scheme at NLO in perturbation the- 
ory at the renormalization scale \i = 1/a [17] , 

(2) B12 are evolved from /x=l/ato/x = 

by using the HQET anomalous dimension matrix, 
known to 2-loop accuracy in perturbation theory 

esusl, 

(3) -Bi, 2 (/i — nib) are then matched onto their 
QCD counterpart, B ia (m b ), in the MS(NDR) 
scheme at NLO [19]. 

The scales chosen to do the matchings are such 
that neither ln(a/x) in the step (1) nor ln(/z/mb) 
in the step (3) correct strongly the matching con- 
stants. The advantage of using a chiral light 
quark action for the step (1) lies in the fact that 
four-fcrmion operators can mix only with a di- 
mension 6 four-fermion operator of the same chi- 
rality. In other words we have not more than 
4 independent renormalization constants in the 
renormalization matrix, because 0\ and 2 can 
mix neither with O3 = h + ^^LS ^~7/iPS, nor with 
4 = h + (l - 7 5 )s hr(l + 7 5 )s. Actually thanks 
to the heavy quark symmetry, those constants are 
not all independent and we have [17j 



8^(M) 
8 2 MS M 



Zu(afi) 

2 (a/j)-Zn(a/j) 





^22 (an) 



4. Results and discussion 



Our results are based on two simulations, with 
the parameters given in Tab. [TJ 
We find 8|^ stat (m 6 ) = 0.922(12) (25), where the 
first error is statistical, the second is systematic 
and contains the error from the estimation of 
a s (l/a) and the finite a effects. From Fig. |3] 
it can be seen that our value is larger than the 
previous static result [5]. This difference is likely 
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1.013(2) 



HYP 


Zl3 


Zu 




no 

yes 


-0.459 
-0.235 


-0.919 
-0.470 


0.763(5) 
0.873(5) 



Table 1 

Parameters of our simulations: amf) and p have 
been chosen following [20121] ; the volume is 16 3 x 
32. 




nd Soni (1998) 



nd Reyes (1999) 



APE (2000) 
DKQCD (2001) 



Table 2 ^ 

Ratios (O 3 ,4) lat /(0i) lat extracted from the simu- 
lation with a Wilson light quark and the improved 
HQET action (left table); comparison of .Z13, 
and i3^ at (TOb) obtained with the Wilson action 
in function of HQET action improvement (right 
table) . 



can take account of the 1/mt effects by interpo- 
lating linearly through Mb s ■ 



B 



MS,A/ B 



(m b ) = B 



MS, stat 



(m b ) 1 + 



C 



(10) 



SPQcdR. (2002) 
JLQCD, NRQCD (2003) 

jlqcd (u „ q . N f = 2). nrqcd (2003) where c _ _q.24(6) GeV [6j. With M Bb = 5.37 

Oraay, static heavy quark (2006) 



0.955(11) x B™ s ' stat (m b ) 



Figure 3. Various lattice values of Bg S (mb) [J- 
[10] ; blue symbols correspond to a computation 
made with a static heavy quark 



due to the use of Neuberger light quark action (no 
subtractions), due to the use of the HYP proce- 
dure, or the combination of both. To answer to 
this question we made a computation with the 
Wilson light quark action. In that case we have 
subtractions in the renormalisation procedure: 



(Oi 



Zn(Oi} lat 



(Q 3 ) lat 
; (Oi) lat 



(Q 4 ) lat 
(Oi) lat 



A correction on Z13 and Z14, coming from a dif- 
ferent definition of the HQET action, induces a 
substantial systematic error on Bb s as illustrated 
on the Tab. [D From Fig. [3] we also notice that 
our value is also somewhat larger than the results 
obtained with the propagating heavy quark. We 



pir nMS,M Bs 

uev, a B ^ — u.auu(nj u Bs 

obtain after the simulation at (3 = 6.0 b^ s < Mb s = 
0.881(15). JLQCD collaboration showed that the 
errors due to quenching are likely to be small 
9 10 . That issue has to be addressed by un- 
quenching the B s — B s mixing amplitude in the 
static limit and by avoiding the subtraction pro- 
cedure as well. The first step would be to make 
a simulation with 2 degenerate Wilson sea light 
quarks and an overlap valence strange quark. 
Eventually with /^ nq = 230 MeV [22], V ts = 
04076 [23], V tb = 0.99912 [24], we obtain 
AM^f = 20.7 ps^ 1 . However this value has to 
be taken very carefully because the uncertainty 
on /™ q is 30%. 

5. Conclusion 

For the first time the bag parameter associated 
to the B s — B s mixing amplitude in the Stan- 
dard Model has been computed on the lattice by 
combining the static limit of HQET and a light 
quark action which preserves the chiral symme- 
try on the lattice. Thus systematic error induced 
by subtractions in the renormalisation procedure 
are absent, since there is no mixing among dimen- 
sion 6 four-fermion operators of different chirality. 
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1 /nib corrections and quenching effects have still 
to be studied carefully. 

REFERENCES 

1. A. Abazov et al [DO Collaboration], 
|hep-ex/0603029| 

2. A. Abulencia et al [CDF Collaboration], 
|hep-ex/0606027| 

3. T. fnami and C.S. Lim, Prog. Theor. Phys. 
65, 297 (1981). 

4. C.W. Bernard, T. Blum and A. Soni, Phys. 
Rev. D 58, 014501 (1998) [hep-lat/9801039] . 

5. V. Gimenez and J. Reyes, Nucl. Phys. B 545, 
576 (1999) [hep-lat/9806023] . 

6. D. Bec irevic et al, Nuc l. Phys. B 618, 241 
(2001) [hep-lat/0 002025 . 

7. L. Lellouch and C. J. Lin, Phys. Rev. D 64, 
094501 (2001) |hep-hep/00lT086| . 

8. D. Becirevic et al, JHEP 0204, 25 (2002) 
and Nucl. Phys. P roc. Suppl. 106, 385 (2002) 
[hep-lat/01 10091] . 

9. S. Aoki et al, Phy s. Rev. D 67, 014506 (2003) 
hep-lat/0208038] . 

10. S. Aok i et al, Phys. R ev. Lett. 91, 212001 

(2003) [hep-ph/0307039] . 

11. K. Hagiwara, S. Narison and D. Nomura, 
Phys. Lett. B 540, 233 (2002) 

12. H. Neuberger, Ph ys. Lett. B 417, 141 (1998) 
hep-lat/9707022] . 

13. E. Eichten and B. Hill, Phys. Lett. B 240, 193 
(1990). 

14. A. Hasenfratz and F. Knechtli, Phys. Rev. D 
64, 034504 (2001) |hep-lat/0103029] . 

15. M. Del ia Morte et al, P hys. Lett. B 581, 93 

(2004) [hep-lat/0307021] . 

16. M. Liischer, Phy s. Lett. B 428, 342 (1998) 
[hep-lat/9802011| . 

17. D. Becirevic and J. Reyes, Nucl. Phys. Proc. 
Suppl. 129, 435 (2004) [hep-lat/03 09131]. 

18. P. Boyle [UKQCD Colla boration], J. Com put. 
Phys. 179, 349 (2002) [hep-lat/9903033] . 

19. D. J. Broadhurst and A. G. Grozin, Phys. 
Rev. D 52, 4082 (1995) [hep-ph/9410240] ; 
V. Gimenez, Nucl. Phys. B 582, 375 (1992); 
X. Ji and M. J. Musolf, Phys. Lett. B 409, 
257 (1991). 

20. L. Giusti, C. Hoelbling and C. Rebbi, 



Phys. Rev. D 64, 114508 (2001); Erratum- 
ibid Phys. Rev D 65, 079903 (2002) 
hep-lat/01080071 . 

21. P. Hernandez, K. Jansen, L. Lellouch and 
H. Witti g, Nucl. Phys. Proc. Suppl. 106, 766 
(2002) [hep-lat/0110"l99l . 

22. S. Hashimoto, Int. J. Mod. Phys. A 20, 5133 
(2005) 

23. M. Bona et al [UTfit group] , [hep^ph/0606167 

24. J. Charles [CKMfitter group], 
hep-ph/ 0606046 



